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1. Introduction 



Precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix is an 
indispensable step to establish the validity range of the standard model, and to search for new 
physics at higher energy scales. Lattice QCD has been making steady progress in this direction. 

One of the difficulties with the charm quark in lattice QCD simulations at a typical cutoff 
' 2 GeV resides in significant cutoff errors due to the charm quark mass. The heavy quark mass 
correction is mga 0{l), and hence it is desirable to control rriQa errors to achieve an accuracy 
of a few percent. The Fermilab action [|l|] and the relativistic heavy quark action ^ have been 
proposed to meet this goal. The relativistic heavy quark action removes the leading cutoff errors of 
0{{mQa)") and the next to leading effects of 0{{mQa)" {aKQCD)) with arbitrary order n, once all 
of the parameters in the heavy quark action are determined nonperturbatively. 

Another source that prevents precise evaluations in lattice QCD is the error associated with chi- 
ral extrapolations in the light quark masses. This problem has been increasingly alleviated through 
progress toward simulations with lighter and lighter dynamical quark masses. The acceleration of 
dynamical lattice QCD simulation using multi-time steps for infrared and ultraviolet modes has 
made it possible to run simulations with light up, down, and strange quark masses around their 
physical values [Q]. In fact, we can proceed one more step and reweight dynamical simulations 
such that dynamical quark masses take exactly the physical values [^. Once the reweighting is 
successfully made, ambiguities associated with chiral extrapolations are completely removed. 

In this paper, we present our work for the charm quark system treated with the relativistic 
heavy quark formalism on the 2+1 dynamical flavor PACS-CS configurations of 32^ x 64 lattice 
generated with the Wilson-clover quark and reweighted to the physical point for up, down and 
strange quark masses. Details of our calculations are presented in Ref. [|6|]. 



2. Setup 



Our calculation is based on a set of A/^/^ = 2 + 1 flavor dynamical lattice QCD configurations 
generated by the PACS-CS Collaboration [||] on a 32^ x 64 lattice using the nonperturbatively 
O (a) -improved Wilson quark action with = 1.715 [0] and the Iwasaki gauge action at /3 = 
1.90. The aggregate of 2000 MD time units were generated at the hopping parameter given by 
(K-o^,K;f) = (0.13778500,0.13660000), and 80 configurations at every 25 MD time units were 
used for measurements. We then reweight those configurations to the physical point given by 
{Kud, Kg) = (0.13779625,0.13663375). The reweighting shifts the masses of 7i and K mesons from 
mn = 152(6) MeV and niK = 509(2) MeV to mj, = 135(6) MeV and mk = 498(2) MeV, with the 
cutoff at the physical point estimated to be a^' =2.194(10) GeV. 

Our relativistic heavy quark action |§] is given by 



Sq = £e,D,,,!2„ 



(2.1) 



[rt - Y4)U,a5,+i,. + {rt + 74)t/l45x,v+4 
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Table 1: Simulation parameters. MD time is the number of trajectories multiplied by the trajectory length. 





Kid 




# conf 


MD time 


1.90 


0.13779625 


0.13663375 


80 


2000 



Table 2: Parameters for the relativistic heavy quark action. 



'Qharm 


V 


rs cb 


ce 


0.10959947 


1.1450511 


1.1881607 1.9849139 


1.7819512 



Cfi L (■^) ^'7 + C£ L ^'4 (■^) 



(2.2) 



where Kq is the hopping parameter for the heavy quark. The parameters r,,rs,CB,CE, and V are 
adjusted as follows. We are allowed to choose = 1, and we employ a one-loop perturbative 
value for [^]. For the clover coefficients cb and ce, we include the non-perturbative contribution 
in the massless limit for three flavor dynamical QCD and calculate the heavy quark mass 
dependent contribution to one-loop order in perturbation theory [^. The parameter v is determined 
non-perturbatively to reproduce the relativistic dispersion relation for the spin-averaged IS states 
of the charmonium. Writing 

E{pf = E{Qf + cls\p\\ (2.3) 

for \p\ = 0, (27r/L),-v/2(27r/L), and demanding the effective speed of light Ceif to be unity, we 
find V = 1.1450511 with which we have Ceff = 1.002(4). It is noted that the cutoff error of 
0{aj{aAQCD)) remains, in addition to C?((aAQco)^), due to the use of one-loop perturbative values 
in part for the parameters of our heavy quark action. 

We tune the heavy quark hopping parameter to reproduce an experimental value of the mass 
for the spin-averaged IS states of the charmonium, given by M(15)'^'^'' = (Mjj^, + 3My/^)/4 = 
3.0678(3) GeV [g]. This leads to JQharm = 0.10959947 for which our lattice QCD measurement 
yields the value M{IS)'"' = 3.067(1)(14) GeV, where the first error is statistical, and the second is 
a systematic from the scale determination. Our parameters for the relativistic heavy quark action 
are summarized in Table Q 

Statistical errors are analyzed by the jackknife method with a bin size of 100 MD time units 
(4 configurations), as in the light quark sector 



3. Charmonium spectrum and charm quark mass 

Our results for the charmonium spectrum on the physical point are summarized in Fig. |l]. 
Within the error of 0.5-1%, the predicted spectrum is in reasonable agreement with experiment. 

Let us consider the 15 states more closely. Since these states are employed to tune the charm 
quark mass, the central issue here is the magnitude of the hyperfine splitting. Our result mjj^^ — 
m-q^ = 0.108(1)(0) GeV, where the first error is statistical and the second error is systematic from 
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Figure 1: Our results for the charmonium mass spectrum(left panel) and the hyperfine splitting of the 
charmonium with different number of flavors(right panel). 



the scale determination, is 7% smaller than the experimental value of 0.117 GeV. In Fig. |l|, we 
compare the present result on A'j = 2 + 1 flavor dynamical configurations with previous attempts 
on Nf = 2 dynamical and quenched configurations using the same heavy quark formalism and the 
Iwasaki gluon action JTOj]. Other results by recent lattice QCD simulations by Fermilab lattice and 
MILC Collaborations |[Tl|], HPQCD and UKQCD Collaborations ^] are also plotted. We observe 
a clear trend that incorporation of dynamical light quark effects improves the agreement. 

We should note that we have not evaluated our systematic errors for the hyperfine splitting, yet. 
The continuum extrapolation needs to be performed. A naive order counting implies that the cutoff 
effects of 0{aj{aAQCD), {<^Aqcd)^) from the relativistic heavy quark action are at a percent level. 
Another aspect is that dynamical charm quark effects and disconnected loop contributions, albeit 
reported to give a shift of only a few MeV [13], are not included in the present work. Additional 
calculations are needed to draw a definite conclusion for the hyperfine splitting of the charmonium 
spectrum. We leave it for a future work. 

Using the axial vector Ward-Takahashi identity, the charm quark mass is obtained. The sys- 
tematic error due to the heavy quark of 0(oi^^(aAgcD), («AgcD)^) will be estimated by data on finer 
lattices in the future. Figure ^ compares our result with a recent Nf = 2+1 lattice QCD estimation 
by the HPQCD Collaboration [|12|]. Another result by ETM Collaboration is also plotted In 
addition to lattice QCD determinations, recent continuum results using the Heavy Quark Expan- 
sions(HQE) [15|, as well as sum rules [16], are shown. All these results are consistent. 



4. Charmed meson and charmed-strange meson spectrum 

We calculate the charmed meson and charmed-strange meson masses which are stable on our 
lattice with the spatial size of L = 2.88(1) fm and a lattice cutoff of a^^ = 2.194(10) GeV. The 
D* and D* meson decay channels are not open in our lattice setup. D*q and Dsi meson masses are 
below the DK threshold [^] but above the D^Tl threshold. Their decays, however, are prohibited by 
the isospin symmetry. On the other hand, Dq and D\ meson masses are not computed since their 
decay channels are open, and therefore a calculation involving Dn contributions is needed. 
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This work(A',=2+l, a-'=2.2 GeV) ^ 
HPQCD(yV/=2+l, a=0) ^ 
ETMCW=2, a=0) ^ 
HQE(BABAR) ^ 
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Figure 2: Comparison of the charm quark mass. We employ Nf =4 running in this plot. 
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Figure 3: Our results for charmed meson masses(left panel) and charmed-strange meson masses(right panel) 
normalized by the experimental values. 



Our results are summarized in Fig. |[ All our values for the heavy-light meson quantities are 
predictions, because the physical charm quark mass has already been fixed with the charmonium 
spectrum. The experimental spectrum are reproduced in 2a level. The potential model predicts 
the D*Q meson mass is above the DK threshold [[l7|], which deviates from the experiment signifi- 
cantly [jl8|]. Our result, however, does not indicate such a large difference from the experimental 



value. A similar result is obtained in other lattice QCD calculations [19, 2^]. The D* 
is below the DK threshold. It should be noticed that our calculation does not cover DK scattering 
states yet. A DK contamination for D*^ and D^i meson masses could be considerably large. Further 
analysis is required to validate our results for D*q and D^i meson spectrum. 

Figure ^ shows our decay constants including the experimental values as well as three 
recent lattice QCD results: HPQCD and UKQCD Collaboration [|T2|], Fermilab lattice and MILC 
group and ETM Collaboration [0]. Our value for fo^ is in accordance with experiment, while 
that for /o is somewhat larger. Comparing four sets of lattice determinations, we observe, both for 
fo and /dj, an agreement between our values and those of the Fermilab group, while there seems 
to be a discrepancy between our values and those by the HPQCD and UKQCD Collaboration and 
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Figure 4: Comparison of pseudoscalar decay constants for the charmed meson (left panel) and charmed- 
strange meson (right panel). 
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Figure 5: Comparison of the CKM matrix elements, \Vcd\ (left panel) and |y„ | (right panel). 



ETM Collaboration, though continuum extrapolation is needed on our part. 

The standard model relates |Vcj| to the leptonic decay width of the D meson r(D — )• Zv) by 

Cj2 / ^2 \ 2 

riD^lv) = ^MmjmD^l-^j \V,,\\ (4.1) 

where Gp is the Fermi coupling constant, and m/ is the lepton mass in the final state. A lattice 
determination of the D meson decay constant fo with the experimental value of r(D — )• Zv) gives 
iKrfl- \Vcs\ can be obtained in the same way. 

We estimate \Vcd\ from our D meson mass and decay constant with the CLEO value of r(D — )• 
Zv) Up to our heavy quark discretization error of 0{aJ{aAQCD),{<^^QCD)^), we obtain \Vcd\- 
Our result of \Vcd\ is about 10% smaller than the PDG value [^, as shown in Fig. ^. Similarly, 
using the CLEO value of r{Ds Zv) [|2l|], we find |V„|, which is consistent with PDG. 
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